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Abstract 



We calculate the special values of the spectral zeta function of the non-commutative har- 
monic oscillator, and give a general formula for them as integrals of certain algebraic functions. 
' This is a generalization of the result by Ichinose-Wakayama (2005), in which the first two spe- 

, cial values are studied. 
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1 Introduction 

> . 

, Let Q be a differential operator defined by 



m 

' The system defined by Q is called the non-commutative harmonic oscillator, which was first intro- 

, duced and investigated by Parmeggiani and Wakayama [9l|T0l[Tl] (see also We always assume 

that a,/? > and a/3 > 1. Under this assumption, Q becomes positive self-adjoint unbounded 
operator on i^(R; C'^) which has only a discrete spectrum. Denote the eigenvalues of Q by 



0<Ai<A2<A3<...(-^oo), 
and define the spectral zeta function of Q by 



This series is absolutely convergent in the region Sfts > 1, and defines a holomorphic function in 
s there. We call this function Cq(s) the spectral zeta function for the non-commutative harmonic 
oscillator Q, which is introduced and studied by Ichinose and Wakayama [T]. The zeta function 
Cq(s) is analytically continued to the whole complex plane as a single- valued meromorphic function 
which is holomorphic except for the simple pole at s = 1. It is remarkable that Cq('S) has 'trivial 
zeros' at s = 0, —2, —4, .... When the two parameters a and /3 are equal, then Cq{s) essentially 
gives the Riemann zeta function ({s) (see Remark [ 
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We are interested in the special values of Cq{^)j that is, the values Cq{s) at s = 2, 3, 4, . . . . The 
first two special values are calculated by Ichinose and Wakayama !_2; as 



Cq(2) 



X (c(2,l/2) 




4:duidu2 



uiui)^ + il-ut){l-ut)/iaP-l)^ 



Cq(3) = 2 



f a + P 
X (c(3,l/2) 



/3\ f 8duidu2du3 



a + Pj y[o,i]3 V(l - uluWa? + (1 - "f)(l - u^ul)/{af3 - 1) 



where C(s, x) — X]^^o("' + 2;)^'* is the Hurwitz zeta function. They also give the contour integral 
expressions of these values using a solution of a certain singly confluent Heun differential equation. 
Later, Ochiai [7] gave an expression of Cq(2) using the complete elliptic integral or the hyperge- 
ometric function (see Remark 15.61) . and the author and Wakayama |S] gave a formula for Cq{3) 
similar to the Ochiai's one. 

In this article, we calculate the special values of the spectral zeta function CqC*^) of the non- 
commutative harmonic oscillator Q for all positive integers n > 1 and express them in terms of 
integrals of certain algebraic functions (Theorem 14. 2p . Roughly, they are of the form 

CQ{n) = 2( " + f ] (c{n,l/2) + Rn{a,(3)], 

and Rn{a, (3) vanishes when a = (3. It is natural to expect that the 'remainder term' i?„(a, P) has 
interesting properties since it reflects the 'non-commutativity' of the operator Q in some sense. 
In fact, for instance, an investigation of R2{a,P) for Cq(2) leads us to a connection with various 
arithmetic objects such as Apery-number-like binomial sums, elliptic curves and modular forms 
(see, e.g., [B]), and we expect that we might have such a kind of connection in the case of i?n(a, P) 
in general. See [3] for the first attempt to study this generalized situation. 



2 Preliminaries 

Following to Ichinose- Wakayama [2], we explain how to calculate the special values of Cq{s). Put 
Bi., V) A- {cos ^i:^/ + sin ^i:^ J 



where 



£;(u,a;,i;) := exp ( -jy(x + y ) + ^—^xy 

l+u 
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and 

^0 13 J ' \p Ij ■■ \1 

Notice that < e < 1 and q > 0. Furthermore, put 

. . . ,x„) := tr iB{xi, X2)B{x2, X3) . . . B{xn,xi)) , 

E{ui, . . . ,U„;Xi, . . .,Xn) ■■= E{ui,Xi,X2)E{u2,X2,X3) . . .E{Un,Xn,Xi), 

F(mi, . . . ,M„) := / E{ui,...,Un;Xi,...,Xn)B{xi,...,Xn)dxi...dXn, 

where tr represents the matrix trace. By Ichinose and Wakayama [5], the integral kernel of the 
operator is given by 



Q-\x,y)^ / A-i/2^(u,x,2/)A-i/2du, 

where 

Kin. x, y) = (1 - e^fl^x. (1 - e^) V^,) exp (^^^ j) • 

Hence, for a positive integer n, we have 
CQ(n) = TrQ-" 

^ y" tr(A-i/2^(Mi,xi,a;2)A-i/2...^-i/2^(-y^^^^^^^)^-i/2j ^^^^ 

^ ^ / F(ui,...,U„) , 



(2.1) 



2 



where dx — dxi . . . dxn, du = dui . . . dun and Tr is the operator trace. Thus, we have only to 
calculate F{ui, . . . ,m„) to get the special values of the spectral zeta function Cq{s). 



3 Lemmas 

3.1 Calculation of the integrand of F{ui, . . . , 
The following lemma is crucial. 

Lemma 3.1. For any positive integer n, it holds that 

B{X1,X2, . . . ,Xn) 



^ ^ ^ [ 0<2k<n ^ I<jl<j2<---<j2fc<n \ r=l 



r 2 
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Proof. Let us put i = ^/—T, ai = a ^, a2 = (3 ^ and tj = qx'j /2 (j = 1, 2, . . . , n). We see that 

B{xi,X2, . . . 



afeiCtfea • • • flfcn COs( - + 

...k„en.2\ m=l ^ 



fei,fe2,...,fc^6{l,2} 



^ ^ ^fcm + l ^mg^(^Tn ^m+l) -|- (^m+1 fcm)g ^(^m ^m+l) 

2j "fci'^fe • • • flfcn n 9 



fei,fe2,...,fc„e{i,2} 

1 

2^ 



A;i,fe2,.-->fen6{l,2} ii,i2,---.'ne{l,-l}™=l 



= - E 

2" 

il,i2,--->'ne{l,-l} fci,fc2,...,fc..e{i,2}"j=i 
= 4 E n(ai^'— '-+a.^^(''^-''"))< 



,i(U-/™-i)t„ 



where we set ko = kn, kn+i = fci, Iq = In, ^n+i = ^i) *o = in and tn+i = ti. Here we notice that 

(ii) |{to e {1, 2, . . . , n} ; Im-i 7^ ^m}| is even (remark that Iq = In), 

(iii) if there exist ji,...,j2/c G {l,2,...,n} such that ji < ••• < j2k, Ij^-i 7^ 'jr fo'^ = 
1,2, . ..,2A; and Im-i = Im for m e {1, 2, . . . ,n} \ {ji, • • • ,i2fe}, then Em=i(^m - lm-i)tm = 

2/,.E'=i(-i)%v- 

Thus it follows that 

B{xi,X2, ■■■,Xn) 

= F 2: (1+ E ir'-a-r\p-'+a-r-'' E eosf2/f;(-l)%.^n 

;e{l,-l} \ 0<2fe<n l<ji<---<j2fc<n V '-=1 / J 

^ \ 0<2fc<n ^ l<ii< - <i2fe<n \'-=l 

This is the desired conclusion. □ 



If we put 



(l-O(l-ut) 








1 4 4 

1— tt;^U2 

(l-ut)(l-u|) 
—U2 





1 44 

1 — UtU^ 



(i-u*)(i-«|) 








1— [t*_j 

.,4 ^ „,4 



(l-<_i)(l-<)/ 
2 



g { (4"' + (t^ - i) + (^a. + ^S'm) t9|j } 
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then 

(3.1) E{ui,. . . . ..,x„) = exp(-a;A„(M)a;') 

and 

(see [U Theorem A.2]). Here we denote by is the matrix unit of size n, and we assume that 

the indices of E",-"^ are understood modulo n, i.e. E^^j = E'^j , ^n+i.j — ^ij^ "^tc. The prime ' 
indicates the matrix transpose. Notice that A„(m) is real symmetric and positive definite for any 
u G (0, 1)". For {ji,j2, ■ ■ ■,j2k} C [n] = {1, 2, . . . ,n}, we also put 

2 k 

s„(ji,...,j2fe):=v^E(-ir4"l- 

r=l 

Since 

2k 

^i~'^y^t = X-nUl, ■ ■ ■,j2k)x' 

and 

we have 

/2k \ 

E{ui, ... ,Un; xi, ... ,Xn) cos gV'(-l)''x|^ 

(3.3) V r=l J 

= 2 (A„(w) + (2'S„(ji, . . . , ja/c)) a;') + 2 ^'^P (^n(w) - '72„(ji, . . . , j2fc)) »;') ■ 

As in [21 Lemma A.l], one prove the 
Lemma 3.2. The determinant 

(3.4) det (A„(m) + qS„(ji, ...,j2k)) 

is even in q. In particular, it follows that this determinant is real-valued for each u G (0, 1)" and 
q>0. □ 

Denote by Cm the cyclic subgroup of the symmetric group of degree m generated by the 
cyclic permutation (1, 2, . . . , m) G &m- By Lemma it follows that 

det (A„(m) + gS„(ji, . . . , jafe)) = det (A„(m) + qE,i{j„(^^, . . . ,ja{2k))) 

for any a G C2k since S„(j^(i), . . .,i„(2k)) = sgn(CT)S„(ji, . . . , ja/c)- 

3.2 LDU decomposition and certain positivity 

Denote by Sym^ the set of n by n complex symmetric matrices such that all principal minors 
are invertible, and by Sym^(R) the set of n by n positive real symmetric matrices. Notice that 
A„(m) G Sym+(IR) for any u G (0, 1)". 
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Lemma 3.3 (LDU decomposition). Let n be a positive integer. For any A e Sym^, there exists a 
lower unitriangular matrix L and a diagonal matrix D such that A = LDL'. Moreover, D is given 
by 

D = diag(rfi,rf2/rfi,d3/d2, . . . ,d„/(i„_i), 
where dk denotes the k-th principal minor determinant of A. 

Proof. Let us prove by induction on n. The assertion is clear if n = 1. Suppose that the assertion 
is true for n — 1. Take A e Sym^ and write 



A = 



Ao a 
a' a 



with Aq G Sym„_j^, a € C" and a G C. By the induction hypothesis, there exist lower unitrian- 
gular matrix Lq and diagonal matrix Dq of size n — 1 such that Aq = LqDqL'q. Put 

where v = {LoDo)^^a and d = a — a'A^^a (notice that (Lo^^o)"^ and Aq^ exist by the induction 
hypothesis) and o € C"~^ represents the zero vector. Then it is straightforward to check that 
A = LDL' . This prove the first assertion of the lemma. The second assertion is obvious by the 
construction of D above. □ 

Lemma 3.4. Let T € Sym^(R) and D be a real diagonal matrix of size n. Denote by dm the 
principal k-minor determinant of T + \/^D. Then it follows that 5i {dm+idm) > for m = 
l,2,...,n-l. 

Proof. Clearly, it is enough to prove the positivity of (dm+idm) with n = m + I. Write T and 
D as 

A a\ „ fU o 



D 



with A G Sym^(]R), a G R™, a G M, m G M and a real diagonal matrix U of size m. Here o G 
is the zero vector. Since T is positive, we must have < a'A~^a < a. Put B = \fA G Sym^( 
X = B-'^UB-'^ G Sym^(R) and h = B'^a. Then we have 



= det 



{A + ^f^U) (A - V^U) a 

a'{A~^/~A_U) a + ^/-iuj 

|det(A + \/^?7)f {a + ^/^u - a' {A + V^Uy^a) 



Since 



it follows that 



or 



= det B^ det(/ + X"^) (a + V^u -h' {l + ^T^X) ^ bj . 
(7 + V^X)~' = (J + X2)-i - V^X{I + xY\ 
^ (b' {L + V^X)~^ bj = b'{L + X^y^b < b'b = a'A-^a < a 



3? 



(a + V^u -b' {1 + V^X) ^ > 0. 



Thus we have ^{dm+idm) > as desired. □ 
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4 Special values 

We recall the well-known fact. 

Lemma 4.1 (Gaussian integral). For any a, 6 G C with 5Ra > 0, it follows that 

exp(^a(x — b)'^)dx — \ —. 

V a 

Here ^/a is chosen as ^^/a > 0. □ 

By Lemma 13.31 ^ G Sym^ is decomposed as A ^ LDL' with a certain lower unitriangular 
matrix L and a diagonal matrix D — diag((ii, d^/di^ . . . , d„/(i„_i), where dk is the A:-th principal 
minor determinant of A. If all entries of D have positive real parts, then it follows from Lemma 
lOthat 



(4.1) / eyiY>{-xAx')dx 

Jb" V det A 

Here x = (xi, .. .,Xn) and da; = da;i • • ■ dxn- 

Now the matrix A„(m) + gS„(j) belongs to Sym,^ for any tt G (0, 1)". Denote by d^ = 
dk{n,u,q,j) the fc-th principal minor determinant of A„(m) +(7S„(j), and put dp = 1- It then 
follows from Lemma [3.41 that di{dk/dk-i) > for k — 1,2, . . . ,n. Consequently, in view of p.3p . 
p.4p . (|4.ip and Lemma [5TT1 we can calculate F{ui, . . . , u„) as 



F(ui,...,'a„) = 2\/^ 



2a(3 



VdetA„(M) oir<„V« + /5y i<„<„t".<,„<„ Vdet (A„(w) + gS„(ji, . . . , j^fe)) 



0<2fc<n ^ ' I<jl<j2< - <j2fc<n 

We also notice that 



[0,1]" 1 - "l"2 ■ • ""n 



= CKl/2) 



for n > 2. From these equations together with (|2.ip and p.2p . we now obtain the 
Theorem 4.2. For eac/i positive integer n > 2, it follows that 

77ere i?„ ^(g) is given by 

2'^dui . . . dUn 

I<jl<i2<---<i2fc<n 



[0,1]" V>V„(M;g;ji,...,j2fe)' 



yV„(M;g; ji, . . . , j2fe) = det (A„(m) + (7^„(ji, . . . , i2k))W{l - ulj) . 



□ 



Remark 4.3. If a = /3, then we have CqC*^) — 2{a^ — 1) "/^C(ri, 1/2), which is a special case of the 
fact that Cq('S) = 2(a^ — 1)~^/^C('5, 1/2) for a = /3. In fact, when a and /3 are equal, we can show 

that Q^V^2-n(^_ij^+ia;2^/ (see [TU]). 
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Remark 4.4. Put uj = y^a/ 13. Consider the limit s — > +0 under the condition that uj is kept 
constant. Then we see from that 



(4.3) 



if n is a positive integer greater than one. Further, for each fixed s > 1, one can prove that 
Cq(s)/£* is analytic with respect to e near the origin and its constant term is (w'* + lli~'^)({s, 1/2). 
The detailed discussion on this result will be given elsewhere [4] . ^ 

We give an expansion of the determinants Wti(w; q; ji, . . . , j2fe) appearing in (|4.2p . For j = 
{jij2, ...,>} C [n] with r > and ji < j2 < ■ ■ ■ < jr, define 



C„(M;j) = n(l 

1=1 

,2, ,2 ,,2^2 



4 4 



i-\-n 



We also define C„(m; 0) = (1 — m^Wj . . . u„) • Here we regard that jr+i = n + ji and u. 
For instance, if n = 9 and j ~ {3, 6, 8}, then 

C9{u;j) = (1 - U3W4u|)(l - M6"7)(l - ulugu^u^). 

Lemma 4.5. for a given subset j — {ji,j2, ■ ■ ■ ijr} C [n] with ji < j2 < • • • < jr, it follows that 

(4.4) Wniu;q;j) = ^(-g^)'^ >V„,rf(M; j) 



d>0 



with 
(4.5) 



W„,,(M;j):= J2 (-l)"^"t^n(w;j(^)). 



SC[2fc] 
|S|=2d 



77ere ||5|| :~ X^sgS ^'^ '^'""^ '^.Z' ^'''^ elements in S and j{S) :— {jsi, ■ ■ ■ ,jsi} if S = {si, . . . , s;} 
with si < ■ ■ ■ < si. 

Proof. Denote by dj the i-th column vector of A„(m). We also denote by {6^}"^^ the standard 
basis of C". By the multilinearity of a determinant, we readily get 

det(A„(M) + gS„(j)) 



2k 



detAr,{u)+J2{^qT (-l)^^+-+^'-det(di,...,e,^^,...,. 



1 (^n 



\<sx<---<s^<2k 



The determinant det(di , . . . , Cj^^ , . . . , e^^^ , . 



, dn) is a product of r tridiagonal determinants 



2a(i) + l,a(i) + l 
^a(i)+2,a(i) + l 



3a(i) + l,a(i)+2 
4(i)+2,o(i)+2 



3o(i) + l,a(i+l)-l 
^o(i)+2,a(i+l)-l 



'^a(i+l)-l,a(i) + l <^a(i+l)-l,a(i)+2 • ■ • C^a(i+l)-l,a(i+l)-l 

where a(i) = j^. , d^j is the (i,j)-entry of A„(m), and the indices are understood modulo n. If 
a(i + 1) = a(i) + 1, then we understand that Di = 1. It is easy to see that 



D, 



1 "a(i)"a(i) + l ' ' ' "a(i+l)-l 



(1 ~ U , — , , ,) ■ ■ ■ (1 — U , , , ) 

Hence we have 

Wn{u;q;j)^ ^ (-1) H^" (^/^q)I^IC„(«; j(5)). 

SC[2fc] 

Since yV„(M; g; j) is real- valued by Lemma IX^ we have the conclusion by taking the real parts. □ 
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5 Examples 

5.1 Wn,diu;j) and Rn,kiq) 

We give several examples of yVn.diu;j). For convenience, we prepare some notation for abbrevia- 
tion. Let us put 

Til / ti \ 

y„(«) := (l-uj. ..ulf, Ut{u) := n 1 - n "'+E.<. ) 

1=1 ^ 3=1 ^ 

for a positive integer n and a sequence t — (ti, . . . , tm) G %n{n), where 

%n{n) := {(ii, . . . ,i„i) G [n]"* ; ti H h i„i = n] . 

For instance, if * = (2, 3, 2, 1) e 7^(8), then 

Ut{ui, . . . = (1 - "l"2)(l - ulu\ul){l - U6W7)(1 - "s)- 

Notice that yV„,o(w; j) = Vn{u) for any ji. 

Example 5.1. For j = C [n] with i < j, we have 

>V„4(m; j) = (-1)^+^C„(m; j) = -(7(r^„_r)(ui, Ui+i, . . . 
where r ^ j — i. This fact immediately implies that Rn^i{q) in ()4.2p is given by 

n-l 



1— ^ /" 

^n.l(9) = ?E / 

2 fr[ J[o 

= E 



2"dM 



0<2r<n 



1]" + 'J'^^nn-r)!") 

n r 2"dM 

1 + S2r,n J[0,1]" a/KI") + (J^ C/(„_r,r) (") 



uf) 



<> 



Example 5.2. For j C [n] with = 2k, it follows in general that 

>V„,fc(M;j) = (-l)^C„(M;j) 
since |1[2A:]|| = fc(2fc- 1) EE A: (mod 2). 
Example 5.3. For j = {ji, j2, js, Ja} C [n] with ji < j2 < h < ji, we have 

WnAu-j) = (-1)1+2C„(m; ji, j2) + (-1)^+^C„(m; Ji, J3) + (-l)i+^C„(M;ii, J4) 

+ (-1)2+3c„(m;j2,J3) + (-1)'+^C„(m;j2,J4) + (-1)3+^C„(m; js, J4) 

= -(1 - 4 • • ■ - 4 • • ■ - 4 • • ■ - 4 • • ■ 

a 4 4 4 4\/i 4 4 4 4\ 

By Example 15.21 we also see that 

yVn,2{u;j) = Cn{u;j) 

= (1 - 4 . . . ui_,)il ~u%... u%_,){l ~u%... - 4 . . . 

Thus we have 



det ( A„ (m) + qS„ ( ji , j2 , is , j4 ) ) n ( 1 

1=1 

= + {q^ + q%l - 4 . . . - U% . . . - 4 . . . - < . . . 
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If wc take (ii, i2, ^3, ^4) G 14(71) such that j^+i — ji + U (mod n) [i — 1,2,3,4; J5 — ji), then it 
foUows that 

du 



[0,1]" ^/det{A 



'[0,1]" V^nH +g^C/(ti+t3,t2+t4)H + +'?'^)C^(tl,t2,t3,t4)(") 

The cychc group C4 of order 4 naturaUy acts on Ti{n) by 

Notice that the integral above is C4-invariant. For a given t = (^1,^2,^3,^4) G ^("■), the number 
of subsets j = {ji, j2, js, ji} in satisfying the condition j'^+i = ji + ti (mod n) is equal to 
n/ \ C4,{t)\, where Ci{t) denotes the stabilizer of t in C4. Consequently, 



T'du 



teTi(n)/C. 

n X ^ 

4 ^ 



|C4(t)| 7(0,1]" V^„(M) + g2J7(i 

+ (g^ + g'')c/(ti,t2,t3,t4)N 



where t — (^1,^2,^3,^4)- Similarly, the result in Example 1 5 . 1 1 can be also rewritten as 

2"du 



teT2{n) 



'[0,1]" V^„(m) +g2(7(ti,t2)(M) 



5.2 Several special values 



<> 



Using Theorem 14.21 and the formulas for i?„^i(q) and Rn,2{<l) given in the previous example, we 
show several examples of the special values of C,q{s). 

Example 5.4 (Ichinose-Wakayama's result). The values Cq(2) and Cq(3) are given by 



<„(3) = 2 ( „ / + ., 1 ( C(3. 1/2) + f 



with 

Aduidu2 f Aduidu2 



i?2,l(<z) 



/[o,i]2 ^wu)Tqnj^:;^ 7[o,i]2 v(i-wM)'+'?2(i-wi)(i-"i)' 

/■ 8duidu2du3 „ /■ 8duidu2du3 

JIO.113 a/V3 M +(7^(7f2.lUM JIO. 



/[0,1]3 ^y3(M)+'Z'C/(2,l)(M) y[0,l]3 v/(l - WlW2"3)' + " ut)(l " U^ut) ' 

This recovers the result obtained by Ichinose and Wakayama [2]. ^ 
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Example 5.5. The values Cq(4) and Cq(5) are given by 



with 



16(itt „ f 16du 



f 32du f 



1]4 ^F4(m) +(72(7(3 7[o i]4 ^F4(m) +g2{7(2,2)(w)' 

32du 



and 



i?4,2(g) 



[0,1]5 VV'sl") +g^C^(44)(^*) "'[0,1]= V^5(w) + q^U(^-i,2){u) 



[o,i]5 ^yVsiu) + (72(7(3^2) N + (9^ + q^)U(^2,i,i,i){u) 



Remark 5.6 (Apery-hke numbers and the Heun differential equation). If we define the numbers 
J2(m) (m > 0) by the expansion 



^2,i(9)=E( i^^)^2(m)<7^ 



2m 



, TO 

m— 

then they satisfy the three-term recurrence relation 

,2 7 /™\ /0^2 Q„, I o\ T /™ T\ I /I/'™ tA2 



Am\J2{m) - (8m^ - 8m + 3) J2(m - 1) + 4(m - 1)^ J2(to - 2) = (to > 1). 
This implies that the generating function W2{t) — X]m=o >^2('Ti)t™ satisfies 

(5.1) |i(i - 1)2^ + (1 _ 3t)(i - 1)| + 1 - ^1 W2{r) = 0, 

which is a singly confluent Heun differential equation 2J . Fortunately, this equation can be reduced 
to the Gaussian hypergeometric defferential equation by a suitable change of variable and solved 
as follows [7]: 



l-t V2 2' 'i-1 

from which we obtain 

i?2,i(g) = 3C(2)2Fi Q,^l;-g2 
Thus we have the following formulas for Cq(2) [H[7]: 

^'^^^^ = (2Va/3(a/-l)) ^ W=T (^) 1(17^4^'^"^ 
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where u{z) — W2(z)/3C(2) is a normalized (unique) holomorphic solution of (|5.ip in \z\ < 1 and 
< r < 1. We also have similar formulas for (q{3) [HE]- Furthermore, if we define J„(m) by the 
expansion of the integral 

i[04]" v/^„H+<7^C/(„-ia)(«) „^o V ™ / 
appearing in _R„.i(g), we also have the recurrence relation 

Am^Jn{m) - (8m^ - 8m + 3) J„(to - 1) + 4(to - l)V„(rn - 2) = Jn-2{m) (to > 1, n > 4) 
and the differential equation 

{t(l - t)^^ + (1 - 3t)(l - 0| + t - ^} .„(t) . (n > 4) 

for Wn{t) — X]m=o (™)^™ ^ ■ However, of course, these results are still insufficient to obtain 
simpler and/or more explicit formula for CqC^*) for n > 4. ^ 
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